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Abstract 

We have investigated the current-induced spin transfer torque of a 
ferromagnet-insulator-ferromagnet tunnel junction by taking the spin-flip 
scatterings into account. It is found that the spin- flip scattering can induce 
an additional spin torque, enhancing the maximum of the spin torque and 
giving rise to an angular shift compared to the case when the spin-flip scat- 
terings are neglected. The effects of the molecular fields of the left and right 
ferromagnets on the spin torque are also studied. It is found that /If> {j^ 
is the spin-transfer torque acting on the right ferromagnet and Ig is the tun- 
neling elcetrical current) does vary with the molecular fields. At two certain 
angles, /le is independent of the molecular field of the right ferromagnet, 
resulting in two crossing points in the curve of /le versus the relevant 

orientation for different molecular fields. 
PACS numbers: 73.40.Gk, 73.40.Rw, 75.70.Cn 
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The spin-polarized transport in multilayer structures exhibits new effects such as the 
giant magnetoresistance [1] (GMR), the spin transfer effect [2], and so on. How to use the 
spin degree of freedom of electrons in ferromagnetic materials to construct new devices is 
at present a focus in the field of spintronics. Spin-polarized electrons flowing from one fer- 
romagnetic layer into another layer in which the molecular field deviates by an angle may 
transfer the angular momentum to the local angular momentum of the ferromagnetic layer, 
thereby exerting a torque on the magnetic moments (see e.g. Refs. [2-9]). This phenomenon 
is usually called the spin transfer effect [2]. The torques in the plane spanned by Si and S2, 
where Si and S2 are spin moments in the left and right ferromagnets, are normally called 
the dynamic nonequilibrium spin torques [9]. Spin transfer motion of Si and S2 within 
their spanned plane is different from the spin precession like dsi/dt — hJsi x S2 out of 
the spanned plane which describes the conventional exchange coupling [12,13]. Therefore, 
the spin transfer effect causes new physical phenomena in magnetic multilayer structures. 
When the current is large enough, it could switch the magnetic states of the local angular 
momentum. Such a current-induced change of the magnetic state has been observed in 
several experiments (see e.g. Refs. [4-7]). As a result, the spin transfer effect may pro- 
vide a mechanism for a current-controlled magnetic memory element. To deal with the 
spin transfer effect, it is useful to introduce the concepts such as the spin current and the 
spin torque to describe the coupling between the conduction electrons and the magnetic 
moments of ferromagnetic materials. These concepts are first proposed by Slonczewski [3] 
based on a quantum-mechanical model for ferromagnet-insulator-ferromagnet (FM-I-FM) 
junctions. Then, the concepts are extended to the structures such as ferromagnet-normal 
metal-ferromagnet (FM-NM-FM) junctions [8,9], ferromagnet-superconductor-ferromagnet 
junctions [10], and trilayer FM-NM-FM which contacts an normal metal lead or a supercon- 
ductor lead [11,12], etc., showing that the investigation on spin torques in magnetic junctions 
has been receiving much attention. 

Many works concerning the spin torque are presented for FM-NM-FM structures so far. 
The work for FM-I-FM structures is still sparse. In particular, when electrons tunnel through 
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the insulator barrier, the spin-flip scattering may occur [14,15]. The spin-flip electrons feel 
a different torque with respect to those non-flip electrons, and could exert an additional 
torque to the ferromagnet. Consequently, this additional torque induced by the spin-flip 
electrons may play a role as the dynamic spin transfer torque. In this paper, we shall use 
the nonequilibrium Green function technique [16] to investigate the spin-flip scattering effect 
on the current-induced spin transfer torque in FM-I-FM tunnel junctions. 

The system under interest is composed of two ferromagnets which are stretched 
to infinite separated by a thin insulator, as illustrated in Fig.l. The molecular field 
in the left ferromagnet is assumed to align along the z axis which is in the junc- 
tion plane, while the orientation of the molecular field in the right ferromagnet, along 
the z' axis, deviates the z axis by an angle Q. The electrons flow along the x 
axis which is perpendicular to the junction plane. The Hamiltonian of the system 
is 

H^Hl + Hr + Ht, (1) 

with 

HL = J24aaiaaka: (2) 
ka 

Hr = - (yM2 cos e)c\^Cq^ - Mssin^cJ^c,^], 

qa 

kqaa' 

where 0^0- and Cfeo- are annihilation operators of electrons with momentum k and spin a 
{— ±1) in the left and right ferromagnets, respectively, e\.^ — ££(k)— e\^— crMi, Mi = ua^a^^ 
M2 — sn^hR^ g is the Lande factor, //s is the Bohr magneton, hL(R) is the molecular field of 
the left (right) ferromagnet, eL(ij)(k) is the single-particle dispersion of the left (right) FM 
electrode, V is the apphed bias voltage, T^^' denotes the spin and momentum dependent 
tunneling amplitude through the insulating barrier. Note that the spin-flip scattering is 
included in Ht when a' = a = —a. It is this term that violates the spin conservation in the 
tunneling process. 

3 



With the system defined above, let us now consider the spin torques exerting on the 
magnetic moments in the right FM electrode of this magnetic tunnel junction. The spin 
torques, namely the time evolution rate of the total spin of the left or the right ferromagnet, 
can be obtained by ^(si,2(i)) = ^{[H, Si,2(i)]). In Refs. [9,11], the spin torques are defined 
by considering the momentum conservation dsi/dt — I{—oo) — 1(0) [2], where I is the spin 
current (whose definition can be found in Refs. [12,8]). Because there are the spin-dependent 
scatterings caused by the local exchange field inside the ferromagnets, the spin current is 
no longer conserved inside the ferromagnets. Whereas the total spin is conserved, the lost 
spin current is transferred to the local magnetic moments, thereby giving rise to a torque 
exerting on the local magnetic moments of the ferromagnets. So the nonconservation of the 
nonequilibrium spin current leads to a current-induced nonequilibrium torque [12]. We may 
consider this issue in another way, i.e. to investigate the evolution rate of the total spin of 
the ferromagnets. In doing so, one must be cautious to identify the spin-torques implied 
from the total Hamiltonian. In this model one may see that the right ferromagnet gains two 
types of torques: one is the equilibrium torque caused by the spin- dependent potential ( i.e. 
the magnetic exchange interaction), and the other is from the electrons tunneling through 
the insulating barrier from the left side. The latter can be obtained from the tunneling term 
Ht, which is nothing but the current-induced spin transfer torque in the plane spanned 
by and Iir. When the apphed bias is absent, the left and right ferromagnets will only 
undergo the spin torque caused by the spin-dependent potential, and the current-induced 
torques would not appear. Therefore, we can define the current-induced spin transfer torque 
as T = ^{[Ht, S2{t)]). Note that a similar expression in a current matrix form is introduced 
in Ref. [10]. The equilibrium torques caused by a spin-dependent potential will not be 
considered here. The total spin of the right ferromagnet is 

^2{t) = ^ E clcU^-'x,y a (R-^x.), (3) 

kixv 



where R 



cos I — sm I 
sin I cos I J 



, (7 is Pauh matrices and XiJ-iu) is spin states. Note that Eq. (3) 



is written in the xyz coordinate frame while the spins S2 are quantized in the x'y'z' frame. 
Prom Si_2 ~ -^6^^1,2 X (si x S2) with Ig the electrical current, we can judge that the direction of 
the spin transfer torque is along x' direction in the x'y'z' coordinate frame. We can further 
write the spin transfer torque in Eq.(3) as S2{t) = | I] {c^^CkwCOsO — ac^^CkaSinO) — 

ku 

S2x'o COS 6 — S2z'fi sin 6, where S2x'o and S2z'q are x'- and ^'-components of the total spins in 
the x'y'z' coordinate frame in which the spins S2 are quantized. Hence, the current-induced 
spin transfer torque can be obtained: 

^ ^Tr.[G<(£)^iTt]+sin^Rei: ' 

kq " kq 



T 



Rx' 



COS 



^ E / ^Tt^[GU^) ^1 Tl + sin ^ Re 1: / ^rr.[G< (e) ^3 Tt], (4) 



where ai- 



A 


u 




A 




°1 


(Ti = 

























are Pauli matrices, T 



Ti T2 



with the elements 



Ti (i = 1,...,4) being the tunneling amplitudes which are for simplicity assumed to be 
independent of A; and q (namely, T^^, T^^ respectively), Tr^ stands for the trace of 
the matrix taking over the spin space, and Gfg{e) is the lesser Green function in spin space 
defined as 



Gll'^ie) G 



kq K^) 



\ 



(5) 



viT,< 

'kq 

with Gl^^''<{s) = / ciie-(*-*')G^f <(i - 1'), and G^;''<(i - 1') = i{cl{t')ak.it)). By using the 
nonequilibrium Green function technique [17], we can get the torque, r^', to the first order 
of the Green function by 



= Try de[f{e + eV) - f{e)\Tra[K{ ai cos 9- sin^)]. 



(6) 



where f{x) is the Permi function, the matrix A is defined as A = T^D£,(£ -|- 



e1/)TRDj^(£)Rt 



D 



Ai A2 
A3 A4 



, and Di(R) 



D 



V 



L(R)T 

D 



whose elements 



L(R)l J 



Lmu) 



[e) = Dn^ii){e ± Mi(2)) are the density of states (DOS) of electrons with spin up 



and down in the left (right) ferromagnet, respectively. After some algebras, we have 

" i / "^^^-^^^^ ~ + eV)](D«t + ^Ri)r[(Pisin^ - Pacos^), (7) 
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.Rx' 



wnere n - z)i^(T^+r|)+Dii(T|+r|)' ^3 - DiT(Tf+T|)+i?i^(T|+Ti') ' and i 1 - i^LTl^i + ^2 J + 
Dli{T^ + T4). It is interesting to note from the above equation that the direction of the 

spin torque is closely related to whether the applied bias is positive or negative, namely, 

it depends strongly on the direction of the electrical current, in agreement of the previous 

observation [2,4-7]. 

To compare the results with those reported in Refs. [2,9], we can also consider the spin- 
torque per current, i.e. t^' /I^ — , where G is the tunnehng conductance [see Eq. 
(12) in Ref. [17]] and V is the applied bias. Then we obtain 

T^' _n Pi sin ^ - P3 cos e 

~ e l + PslAcos^ + Pgsin^)' 
where P2 = is the polarization of the right ferromagnet, and the energy is taken at 

the Fermi level. 

To gain deeper insight into the effect of the spin-flip scatterings on the spin tranfer 
torque, we ought to invoke numerical calculations. Before presenting the calculated results, 
we shall presume a parabolic dispersion for band electrons based on which the DOS of 
conduction electrons are calculated. The Fermi energy and the molecular field will be taken 
as Ef = 1.295 eV and |hi| = |h2| = 0.90 eV, which are given in Ref. [18] for Fe. In addition, 
we may for convenience introduce two parameters 71 = T2/T1 and 72 = T3/T1, and assume 

Now let us first look at the case without the spin-flip scatterings. The leading contribu- 
tion to the torque t^^' comes from the first nonvanishing term. When the spin-flip scatterings 
are absent, i.e. T2 ^ ^ 0, we find that = f (^RT + Dri){Dl^ + D^jT^Pi sin ^ 

with Pi = {Dli — Dli)/{Dli + -Dlj,), which is consistent with that obtained in Refs. [2,3]. 
This result shows clearly that t^'"^' vanishes when the relative alignment of magnetizations 
of the two ferromagnets is parallel {9 — 0) or antiparallel {9 — tt). It is similar to that 
for a FM-NM-FM trilayer system discussed in Ref. [9] , although the transport mechanisms 
are different. The present result can be easily understood because the spin-polarized elec- 
trons along the z or —z axis cannot feel the spin transfer torque owing to the property of 
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si,2 X (si X 52)- In Fig. 2, we show the 9 dependence of the spin torque in the absence 

Rx' 

of spin-flip scatterings for different polarizations. One may see that is a nonmonotonic 

Rx' 

function of 9, and shows minima and maxima at certain relative alignments. versus 9 is 
inversion-symmetrical to the axis of ^ = 0. It is evident that the larger the polarization, the 
stronger the spin transfer torques. This observation is in good agreement with the finding 
in Ref. [2], though the latter is obtained on a basis of a quite different method. 

Rx' 

The ^-dependence of for 71 = 72 is shown in Fig. 3 (a). It can be seen that the spin- 
flip scatterings can lead to a nonvanishing spin torque at ^ = or tt, i.e. ^'^^ = =pe7r(Di^| -|- 
^Rl)Ti{'yiDL-^ + j2DLi), being different from the case when the spin-flip scatterings are 
neglected. However, in the present case with 71 = 72 = 0.05, at a particular angle, e.g. 
9 — 14.04°, T^^' becomes zero. This suggests that the spin-flip scattering can lead to 
an angular shift to the spin torque. In other words, the spin-flip scattering gives rise to 
an additional torque which we may call the spin-flip induced spin torque henceforth. In 
addition, if 71 = 72 = 7, one may find that gy is propotional to 7. When 71 7^ 72, 
which means that the spin-flip scatterings from the spin up band to the spin down band are 
different from those from down to up, 71 and 72 give different effects on the ^-dependence of 
spin torques as shown in Fig. 3 (b). It is observed that there are angular shifts for different 
7i and 72, showing that the effects of 71 and 72 are various. For instance, the curve for 
7i = 0.2 and 72 = 0.1 moves to the right-hand side in comparison to the curve of 71 = 0.1 
and 72 = 0.2. It appears that a larger spin-flip scattering from spin down to spin up has a 
larger effect on the spin torque. 

The effect of the molecular fields of the ferromagnets on the spin torques is also investi- 
gated. The ^-dependences of r^' are shown in Fig. 4 for different parameter a — Ih^l / Ih^l. 
One may find that t^' /le is also a nonmonotinoc function of 9, and shows peaks at certain 
values of 9 for different a's. It is interesting to note that the two crossing points at 9i — 37° 
and 92 = 127° under different molecular fields are observed. When ^ < ^1, a larger a leads 
to a smaller magnitude of the spin-torque (which are negative); when 9i < 9 < 92, the larger 



7 



a, the smaller the spin-torques; while 6 > 62, at, larger a leads to a larger magnitude of the 
spin-torque, showing that |h^| and |hi| have different effects on the spin torques. The two 
crossing points appeared in the curve T^^'//e versus 9 for different molecular fields can be 
understood in the following way. Prom Eq. (8), one can get a relation 

h ^ 1 + P2^J p^ + Pi cos{e- Of) 

where the angular shift ^/ is defined by tan^j = ^ [17]. 9i at which the first crossing point 

Rx' 

appears is nothing but 9f (i.e. 9i — 9f — 37°). At 9i = 9f, = 0, which is independent of 
the molecular field of the right ferromagnet. The second crossing point appears at 92— 9 f — | 
(i.e. 92 — 127°). At this angle 92, ^^^^ = f + Pi, which depends only on the parameters 
of the left ferromagnet but is independent of the molecular field of the right ferromagnet. 

Rx' 

Therefore, at these two particular alignments, ^^j— gives the same value for different a's, 
thereby leading to the two crossing points. When the spin-flip scattering disappears, the 
two crossing points occur at ^ = and |. So we may find that the spin- flip scattering may 
lead to an additional spin torque on the magnetic moments of the right ferromagnet. The 
present observation may be readily understood, because the spin-flip tunneling of electrons 
from the left ferromagnet into the right ferromagnet would give rise to a change of the DOS 
of conduction electrons as well as the effective polarization factor in the right ferromagnet, 
leading to an additional torque exerting on the magnetic moments, thereby the property of 
an angular shift of the spin torque was observed. 

In summary, we have investigated the current-induced spin torque in FM-l-FM tunnel 
junctions with inclusion of spin-flip scattering by using the nonequilibrium Green function 
method. In the absence of the spin-flip scattering, our results are consistent with the previous 
results found in Refs. [2,3,9]. When the spin-flip scattering, the factor that could exist in 
realistic spin-based electronic devices, is considered, we have found that an additional spin 
torque is induced. It is found that the spin-flip scattering can enhance the maximum of 
the current-induced spin transfer torque, giving rise to an angular shift. The effects of the 
molecular fields of the left and right ferromagnets on the spin torques are also studied. 
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It can be observed that the spin-torques per unit tunnehng current exerting on the right 
ferromagnet are independent of the molecular field of the right ferromagnetic lead at 9 — 9f 
and — 9f = |. The present study shows that the spin- flip scatterings during the tunneling 
process in magnetic tunnel junctions have indeed remarkable influences on the dynamics 
of the magnetic moments. Since the spin-transfer torque induced by the applied current 
can switch the magnetic domains between different orientations [6], people can invoke this 
property to fabricate a current-controlled magnetic memory element. As the previous studies 
on the spin torque usually ignore the effect of the spin-flip scatterings which could not be 
avoided in practice, our investigation might offer a supplement to the previous studies, 
namely, when people design a device based on a mechanism of the spin transfer torque, one 
should take the additional torque induced by the spin-flip scatterings into account, which 
could complement with the experiments. 
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FIGURE CAPTIONS 

Fig. 1 A schematic illustration of the spin-transfer torque with spin-flip scatterings in a 
FIF tunnel junction. 

Rx' 

Fig. 2 Spin transfer torque versus angle 9 in the absence of the spin-flip scatterings 
(i.e.7i = 72 = 0). Equal polarization factors of the magnets are assumed (Pi = P2). Torque 
per unit current is measured in unit of h/e. 

Fig. 3 The 6 dependence of for (a) 71 = 72 and (b) 71 7^ 72, where the effective 
masses of the left and right ferromagnets are taken as unity, the molecular fields are assumed 
to be 0.9 eV, the Fermi energy is taken as 1.295 eV, and Ti — T4 — 0.01 eV. Torque per 
unit current is measured in unit oih/e. 

Fig. 4 r^^'/-^e versus angle for different molecular fields. Here we take 71 = 72 = 0.15, 
|hi| = 0.9 eV, and the other parameters are taken the same as those in Fig. 3. Torque per 
unit current is measured in unit oi h/e. 
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